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Mineral Dust Module in GEOS-Chem 


• Dust mobilization (Fairlie et al., 2006) 

•Transport (TPCORE/FvDAS) (Lin and Rood, 1996) 

•Dry deposition - gravitational settling, turburlent transfer, 
impaction, interception, Vd=Vd(D,r,sfc) 

(Zhang et al. 2001) 

•Wet deposition - scavenging in convective updrafts; 
rainout and washout from large-scale precipitation and 
convective anvils (Liu et al., 2001). 

•Size bins: 0.1 -1.0, 1.0-1. 8, 1. 8-3.0, 3. 0-6.0 pm radius 
(following Ginoux et al., 2001) 



Mobilization 


Vertical dust flux, Fd, proportional to the horizontal 
saltation flux, Qs: Fd ~ S Qs 

• Qs ~ U 10 2 (U 10 - U* t ) (GOCART) 

• Qs ~ U* 3 (1 -U7U*) (1 +lT t /U*) 2 (DEAD) 

U* t - threshold friction velocity (particle size, density; air 
density, viscosity). U* t modulated by surface moisture. 

S - source function (defines potential dust source 
regions, and comprises surface factors, e.g. 
vegetation, snow cover, and an efficiency factor e.g. 
topographic anomaly) 


Emission, average load, wet and dry 

depositions for 2001 
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Observed and simulated mineral dust during INTEX-B: transpacific 
transport and evidence of acidic uptake on dust surfaces 


We use measurements of aerosol ion composition made from the 
SAGA instrument (Dibb et al., UNH) on board the DC8 aircraft 
during the 2006 INTEX-B airborne campaign to identify mineral dust 
signatures. 

Use Na+ and Ca(2+) to distinguish sea salt and mineral components 
of the aerosol distribution, respectively. Positive correlations of non- 
sea-salt sulfate and nitrate with calcium indicate that the dusty air 
stream is associated with polluted air masses. 

Look for evidence for interaction of dust with acidic components: 

Coating of dust with sulfate or nitrate (i) shifts sulfate or nitrate 
towards larger sizes, (ii) favors uptake of water on dust surfaces, 
promoting dust particles as cloud condensation nuclei, (iii) promotes 
further uptake of S0 2 and of N 2 0 5 , (iv) impacts N0 X /HN0 3 
partitioning, and (v) affects atmospheric lifetimes for both aerosol 
and gas components. 



NTEX-B (SAGA): 20060417 - 20060515, Pacific 
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Scatter plots of Ca(2+), Mg(2+), and K+ vs. Na+; Mg(2+) and K+ vs. Ca(2+), from the 
SAGA bulk aerosol measurements for the Pacific deployment of INTEX-B. Blue symbols 
indicate observations made below 1 km altitude, showing distinct sea salt signatures 
(blue lines). Above 1 km, observations show different but equally distinct signatures: 
Ca/Na - 2.5; K/Na - 1/4; Mg/Ca - 3/10. For Mg(2+) and Ca(2+), r 2 = 90% above 1 km 
with a slope of approximately 0.3. This population with Ca(2+) reaching in excess of 
1200 ppt is characteristic of mineral dust. 


20060515, Pacific 


NTEX-B (SAGA): 20060417 
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Scatter plots of S04=, N03-, and HN03 vs. Ca(2+); S04 = and N03 vs. NH4+ from 
SAGA bulk aerosol data for the Pacific deployment of I NTEX-B. [data limited to z > 1 


> S04= is positively correlated with Ca(2+) => the dust is often associated with pollution. 

> For S04= vs. NH4+ , r 2 = 71%, slope = 0.44, => S04 primarily in the form of (NH4) 2 S04. 

> N03- vs. NH4+ shows little evidence for (NH4)N03. N03- positively correlated with Ca(2+) 
(r 2 = 70%, slope = 1/3), while HN03 shows an inverse relationship with Ca(2+). 

> Taken together, this suggests that much of the nitrate is in coarse mode, i.e. on the dust. 

> The slope of N03- vs Ca(2+) is well below 2, leaving excess Ca(2+). 


GEOS-Chem 3-d model results: 


Results shown from a full chemistry and a separate dust-only simulation for the 
INTEX-B period conducted at a horizontal resolution of 2 deg. by 2.5 deg. 

The model is driven by GEOS-4 assimilated meteorological fields. 

The simulation represents OC, EC, S04-N03-NH4, sea salt, and dust aerosols. The 
latest version (v7.03.06) includes updated treatment of aerosols: (i) thermodynamical 
equilibrium is now computed using ISOROPIA (Nenes et al.); (ii) coarse mode S0 4 and 
N0 3 on sea salt aerosols (S0 4 s, NITs) (Alexander et al.); (iii) hygroscopic growth of 
aerosols is now accounted for in dry deposition; the emission, transport, wet and dry 
deposition of mineral dust (Fairlie et al.). Full chemistry results are provided courtesy of 
Lin Zhang, Rokjin Park, Bob Yantosca. 

To compare modeled dust with SAGA data, we computed Ca(2+) from the dust 
component, following Song and Carmichael (2001) in assuming Ca(2+) comprises 6.8% 
of (Asian) dust by mass. 

NB: Comparisons of simulated and observed Ca(2+) indicate that simulated dust 
concentrations are 4 times too large. The results shown here show simulated dust 
(Ca2+) scaled by a factor of 1/4. 



Simulated April/May 2006 dust emissions, average load, dry and wet depositions 

[145Tg; 14.5 Tg; 53 Tg; 94 Tg, respectively 
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Dust AOD and surface concentration averaged for the Pacific deployment (20060417- 
0515). Contours of average mean sea level pressure are also shown. The map shows 
the mean transpacific transport of mineral dust from Asia dust storm activity to North 
America during April-May 2006. SAGA data points along the DC8 flight tracks (black 
dots) were well positioned to sample the transpacific dust transport. 
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Dust mixing ratio at 40°N averaged for the INTEX-B Pacific deployment (20060417- 
0515). The section highlights the principal source regions and the mean transport across 
the Pacific during April-May 2006. In addition to dust emissions from the Gobi and 
Taklimakan deserts, the model indicates upstream contributions to transpacific transport. 



NTEX-B Pacific 20060417-0515 
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Latitude-altitude cross-sections of mean simulated dust mixing ratio and eastward dust 
flux at 180°W during the Pacific deployment (20060417-0515). Isopleths of potential 
temperature and zonal wind speed are also shown. The sections show the bulk of the 
dust between 30 and 60°N, with the largest eastward flux around 40°N. SAGA data 
points along the DC8 flight tracks are also shown (black dots). 


INTEX-B (SAGA) 
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20060417 - 20060515, Pacific 
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Altitude distribution of >AGA ion data vs. model data. 

> Medians of simulated S04 and NH4 generally fall between the quartiles of the 
observations throughout the troposphere. 

> Simulated N03 is biased low in the low to mid troposphere, whereas HN03 is 
biased high throughout the troposphere. 

> Modeled Ca(2+) shows the dust peak between 1 and 3 km whereas SAGA 
indicates the peak between 4 and 5 km. 


NTEX — B (GEOS — Chem): 20060417 
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Scatter plots from the model corresponding to those shown earlier for the SAGA data. 

> Like SAGA, the model shows a positive association between S04= and Ca(2+). 

> S04= vs. NH4+ shows a slope of about 0.5 indicating (NH4) 2 S04. 

> Some indication of excess NH4+ and evidence for (NH4)N03, in contrast with SAGA. 

> Positive association of Ca(2+) with HN03, in contrast with SAGA. 

> Lacks of the positive association between N03- and Ca(2+) that was found in obs. 

> Differences between the simulated and observed N03 partitioning indicates the 
importance of HN03 uptake on dust. 


Simple representation of uptake of HN03 on dust: 


Ca(C03)(s) + 2HN03(g) Ca(N03)2(aq) + H2C03(aq) 

-> Ca(N03)2(aq) + H 2 O + C02(g) 

d [HN03] /dt = -2K [Ca] a [HNO 3 ] 
d [N03] /dt = 2K [Ca] a [HNO 3 ] 

d [Ca] a / dt = -K [Ca] a [HN03] where [Ca]a is the available Ca(2+); 
Keeping [Ca]a constant yields 

[HN03] = [HN03] 0 exp(-y) 

[N03] = [N03] 0 + [HN03] 0 ( 1 — exp(-y) ) where y = 2K [Ca]a t ; 

[HN03](t+5t) = [HN03](t) -2K [Ca]a [HN03](t) 8t 
[N03] (t+8t) = [N03] (t) + 2K [Ca]a[HN03](t) 8t 
[Ca] a (t+St) = [Ca]a (t)-K[Ca] a [HN03](t) 8t 


Let y = 2K [Ca]a t = 0.0006 



Impact of simple representation of HN03 uptake on dust 
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Impact (i) reduces HN03 at 
high Ca(2+) values 


(ii) raises N03- at 
high Ca(2+) values 


(iii) raises N03- vs. NH4+ 
Closer to observations 


NTEX-B (SAGA): 

N03 vs Altitude 
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Simple scheme: 

> raises N03- mixing ratios in the mid-upper trop. comparable with observations; 
larger between 1 and 4 km coincident with peak dust. 


> We do have a total N03- (= [N03-] + [HN03]) high bias in the model in the trop 


Conclusions 


The SAGA observations indicate sulfate to be primarily in the form of 
(NH 4 ) 2 S0 4 . A positive relationship between Ca(2+) and N0 3 -, together with 
an inverse relationship between HNO s and Ca(2+), and little evidence for 
NH 4 N0 3 , suggest that N0 3 - is associated with mineral dust surfaces. 

3-d model simulations conducted with GEOS-Chem indicate that 
transpacific transport was principally responsible for the dust observed from 
the aircraft over the Pacific. The model shows a very different partition with 
respect to N0 3 - : elevated HN0 3 and substantial evidence for NH 4 N0 3 . 

Simple representation of HN0 3 uptake raises N03- mixing ratios in the 
lower troposphere coincident with the simulated dust peak, and improves 
simulated ion correlations. 


These results indicate that the uptake of HN0 3 on dust surfaces, missing 
from this model formulation, is an important process in accounting for the 
partitioning of N0 3 -in the INTEX-B data. 

Stick around for Fabien’s presentation 



Questions 



[Typically for seawater: Na+ = 0.469; 
Mg(2+) = 0.0528; S04(=) = 0.0283; 
Ca(2+) = 0.0183 (mol/kg)]. 

Dust 

Ca:Na ~ 5:2; K:Na - 1:4; Mg:Ca ~ 3:10 

[Na+ is also a component of dust, Na20, 
typically 2% by weight]. 



